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Abstract. We calculate the evaporation of a cool accretion disk around a black hole due to the ion-bombardment 
by an ion supported accretion flow (here ISAF, or optically thin ADAF). As first suggested by Spruit & Deufel 
( 2002), this evaporation takes place in two stages: ion bombardment of the cool disk (Shakura-Sunyaev disk: SSD) 
produces an intermediate-temperature layer on top of the disk ('warm layer') which constitutes an independent 
accretion flow on both sides of the SSD. As this warm material accretes inward of the inner radius of the SSD, 
it becomes thermally unstable by lack of cooling of photons, and evaporates into the ISAF, thereby feeding the 
latter. Angular momentum conservation forces a certain fraction of the ISAF material to move outward, where it 
can bombard the SSD with its hot ions. The flow geometry is derived by computing stationary solutions of the 
continuity- and angular momentum equations for the three components (ISAF, warm flow and SSD). The overall 
radiative output is dominated by hard X-rays. They are produced mostly from the warm component, rather than 
the ISAF. The expected time dependence and stability of the flow, not computed here, is discussed briefly. 
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1. Introduction 

The X-ray spectra of galactic black holes occur in two 
different 'states': a soft state dominated by a soft multi- 
color blackbody component peaking around 1 keV, and a 
hard state in which the energy output peaks at 100-200 
keV. The soft state can be explained well with a stan- 
dard optically thick geometrically thin accretion disk of 
the kind of Shakura & Sunyaev l|l973jh But what pro- 
duces the emission in the hard state is still under debate. 
Some argue that it represents non -therm al synchrotron 
radiation from a jet (Markoff et al l200ll) . Others invoke 
magnetic activity on the surface of the di sk, cr eating a 
'coro na' on top of the d isk (G aleev et al. Il979t Haardt 
et al. Il994t Stern et al. Il995l Merloni & Fabian 1200 llh 
A commonly cited explanation for the hard X-ray spec- 
tra of black holes is the comptonization of soft photons 
in of a hot optically thin virialized plasma flow. The con- 
cept of such a hot 'ion torus' around black holes has al- 
ready been propo sed a long time ago (Shapiro, Lightman 
& Eardlev 11970: Ichimaru Il977r. Rees et al. Il982|h but 
optically thin plasmas in the 100 keV range are noto- 
riously unstab le (Li ghtma n & Eardle v Il974t Sunyaev & 
Shakura[l97l Piran ll978t Chen et al. ll995lh In the 90s it 
was recognized that the advection of thermal energy in an 
inefficiently radiating hot optically thin plasma fl ow ca n 
have a stabilizing effect on the flow (Narayan & Yi ll995|h 
Instead of converting potential energy into radiation, such 



an "Advection Dominated Accretion Flow" (ADAF, here 
called "Ion Supported Accretion Flow" : ISAF) stores the 
released potential energy into thermal energy of the ions 
and simply advects it down through the event horizon. A 
crucial element in the theory of such ISAFs is the ther- 
mal decoupling of ions and electrons. If the temperature 
of the electrons would be equal to the (almost virial) tem- 
perature of the ions, the flow would cool dramatically and 
collapse. The weakness of the Coulomb coupling between 
the eletrons and the ions allows the electrons to have a 
much cooler temperature (of the order of 2 x 10 9 K) than 
the ions (which can have temperatures reaching 10 11 K), 
and therefore to radiate only very limitedly, allowing the 
thermal energy to be advected into the black hole and the 
hot flow to be stable. Despite the relative inefficiency of 
ISAFs, they do emit at least some Bremsstrahlung and 
thermal synchrotron radiation, and they comptonize soft 
photons to high energies, which could explain the hard 
X-ray components observed in the spectra of black holes. 

While this scenario seems appealing, there are a num- 
ber of open issues. One of these is the existence of 'in- 
termediate' states in which both components are present 
simultaneously, suggesting that the hot ISAF and the cool 
SSD somehow coexist. Based on simple mean-free-path ar- 
guments one can show that if these flows exist at the same 
distance from the black hole, the ions of the ISAF (with 
energies of 1 - 100 MeV) are quickly lost by penetrating 
into the SSD, causing the ISAF to condense. Another issue 
is that ISAFs can only exist up to a certain distance from 
the black hole, because the two-temperature plasma state 
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ceases to ex ist at low tem peratu res, (e.g. Igumenshchev 
et al. Il998l Czerny et al. |200(A Narayan & Yi Il995|) . 
Therefore their existence means that they must have been 
produced in situ, probably by some evaporation mech- 
anism acting on the cool Shakura-Sunyaev disk (here- 
after 'SSD'): the ISAF being the 'vapor' from the SSD. 
A well-studied evaporation mechanism is 'coronal' evapo- 
ration (which determines the chromosphere-corona transi- 
tion in the solar atmosphere). It is governed by the conduc- 
tion of heat by the electrons (Meyer & Meyer-Hofmeister 
Il994l) . This has been shown to work well for models of 
Cataclysmic Variables. In the two-temperature regime, 
this process fails because the energy is not carried by the 
electrons any more. For black hole accretion the mecha- 
nism therefore only works well at large radii (> 10 3 i? s ). 
Yet there are strong indications that some evaporation 
mechanism acts on the SSD at much smaller radii: from 
analysis of X-ray spectra of the soft component it has been 
inferred that the inner radius of the SSD regularly moves 
outward from the last stable orbit to distances of a few 
tens of Schwarzschild radii, in particular during the transi- 
tion fr om a 'high-soft' state t o a 'lo w-hard' state (Zd ziarski 
et al. 119991: Gilfanpv et al. l2000l Churazov et al. l200ll 
Revnivtsev et al l200l|) . 

The simultaneous presence of a hot comptonizing 
plasma close to the black hole and an SSD with varying 
inner radius therefore poses a problem: what is the mech- 
anism that evaporates the SSD at close distances from 
the event horizon, thereby prod ucing the hot ISAF? Some 
time ago Spruit & Deufel (2002) proposed that a possible 
mechanism could be in fact the injection of ISAF ions into 
the SSD. Their calculations show that such an ion bom- 
bardment produces an intermediate-temperature (about 
60 keV) layer of optical depth t ~ 1 on t op of the SSD 
surface (Deufel, Dullemond & Spruit 12002^ . 

The temperature of the warm layer is kept stably in 
this intermediate range (which would be thermally un- 
stable for an isolated optically thin plasma) by the bal- 
ance between heating and Compton cooling by the soft 
photon flux from the cool disk. This is similar to the en- 
ergy balance in the Haardt-Maraschi model (1994), but 
the present model is more predictive, since both the tem- 
perature and the optical depth of the layer are determined 
by the physics of ion penetration. 

Such a 'warm layer', being of much lower density and 
much higher temperature than the SSD (but still much 
lower than the virial temperature) , transports its own an- 
gular momentum and forms an accretion flow on top of 
the disk. Now suppose the inner radius of the SSD (i? ss d) 
is significantly far outside the last stable orbit. The ra- 
dial flow velocity in the warm layer is orders of magnitude 
higher than the flow velocity in the SSD, and therefore 
this warm material quickly accretes inward of i? S sd- In 
this region there is no SSD underneath the warm layer 
to provide soft photons for compton cooling, and there- 
fore the temperature of the warm accretion flow inwar d 
of i? ss d increases (to roughly 300 keV, Deufel et al. l2002l) . 
Due to the reduced soft photon flux, this region is now 
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Fig. 1. The geometry of the three-layered accretion flow, 
in which evaporation of the cool disk takes place in two 
stages. 



also thermally unstable, and starts evaporating into the 
ISAF (Spruit & Deufel 2002) 

This process produces the mass flux that sustains the 
ISAF. Part of the ISAF material accretes inward, and part 
moves outward due to angular mo mentu m conser vation 
(as in any viscous disk, e.g. Pringle 1981; Liu et al. 1997; 
Turolla & Dullemond l2000j) . The outward moving ISAF 
material spreads over the SSD. Over a short distance it 
dumps its ions onto the SSD, thus producing the warm 
layer and closing the cycle. This proposed mechanism is 
therefore a two-stage evaporation process: first from the 
SSD into the warm layer, and then from the warm layer 
(at smaller radii) into the ISAF (see Fig.^for a sketch). 

In the paper of Spruit & Deufel 'zero dimensional' 
estimates were given for the three-component flow ge- 
ometry, without actual computation of the flow geome- 
try. This is the topic of the present paper. We solve the 
combined continuity- and angular-momentum conserva- 
tion equations for all three components (ISAF, warm flow 
and SSD) simultaneously with the appropriate mass ex- 
change terms between them. We derive stationary solu- 
tions to these combined equations and study their prop- 
erties. 



2. Basic model 

We treat all three components (ISAF, warm flow and SSD) 
in a vertically integrated way, but we allow the compo- 
nents to coexist with each other at every radius. Therefore, 
each component is represented by a surface density func- 
tion S x (-R, t) where R is distance from the center, t is time 
and x stands for either "ISAF" , "warm flow" or "SSD" . 
Both the warm flow and the SSD may be truncated at 
an inner radius (larger than the last stable orbit). These 
radii, called i? W arm for the warm layer and R ss d for the 
SSD may vary with time. In the stationary state they are 
eigenvalues of the problem and depend on the global ac- 
cretion rate. The temperature of the ISAF is assumed to 
be a constant fraction times the local virial temperature: 



T isaf = 1.26 



VI + 0.56 a 2 - 1 GM ^m p 



R 



(1) 
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which can be derived from self-similar solutions of ADAFs 
using 7 = 1.5 (Narayan & Yi ll994|) . The temperature of 
the SSD is calculated self-consistently assuming Kramer's 
opacity: 



T ssd (R) = 1.17 x 10 8 aS 



1 M/Md 



1/7 



(2) 



The temperature of the warm layer is assumed to be 300 
keV for R < R ssd and 60 keV for R > R ssd . These val- 
ues are e stimates based on the results from Deufel et 
al. (|002), who calculated the detailed vertical structure 
of the warm flow. 



2.1. Evaporation of warm layer into ISAF 



Spruit & Deufel ( 2002) derived the evaporation rate of the 
warm layer inward of i? ss d into the ISAF. They find the 
following expression (valid for R < R ss d)- 
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(where kT e w 300 keV). The evaporation rate is there- 
fore constant between the inner radius of the warm layer 
and the inner radius of the cool disk. The inner radius 
of the warm layer i? war m may be located at the last sta- 
ble orbit at R = 6GM/c 2 , but it can also be larger if 
the evaporation process was so effective that it evapo- 
rates the layer before it has reached the black hole. At 
this stage, we assume the latter, i.e. that the warm layer 
has i?warm > QGM/c 2 . The consistency of this assumption 
will be part of the results obtained. 

2.2. Evaporation of SSD into warm layer 

The evaporation from the SSD into the warm layer for 
R > i? S sd is a more subtle matter. This evaporation is 
powered by the ion bombardment of the SSD by the ISAF. 
The flux of ions from the ISAF onto the SSD (and through 
the warm layer) is: 



isaf — >ssd 
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which is defined as the sum of the fluxes on both sides 
of the disk. The parameter £ is a dimensionless efficiency 
parameter of order unity. If £ = 1, the ISAF loses its 
matter at the fastest possible rate: any proton that reaches 
the disk surface is lost. Effects such as tangled magnetic 
fields or proton-proton collisions may reduce the loss rate 
of matter to the disk, which is simulated by taking £ < 1. 

The protons that enter the SSD, leave their energy 
in the disk's upper layers. A fraction rj of this energy is 
radiated away by Bremsstrahlung and possibly other pro- 
cesses. It is not straightforward to compute this value T], 
so we keep it a parameter. The remainder of this energy 



can be used for the heating of disk matter up to the tem- 
perature of the warm layer. This creates an evaporation 

rate 'ipssd — »warm of. 

^isaf -Mvarm 



V'ssd- 



(1 - ry)^isaf^ssd: 



(5) 



From the static solutions of Deufel et al. (2002) we 
know that, as S approaches a value close to unity (let us 
define this value as So), the evaporation will cease and 
a static solution sets in. In this static solution the en- 
ergy input by ion bombardment is entirely compensated 
by radiative losses (i.e. radiative losses in excess of the 
•q fraction already taken into account in the above evap- 
oration rate formula). We simulate this saturation effect 
by multiplying ^ssd^ssd by a factor (1 - /S ). The 

mass gain of the warm layer is then: 

^ssd^warm — 

Cd —►warm I v 1 

- ■ (6) 

Here, and in the remainder of this paper, we take So = 6 
g/cm 2 (accounting for two sides of the Shakura-Sunyaev 
disk, each having a layer of maximum surface density of 3 
g/cm 2 .). 

2.3. Equations for all three components 

All three layers (disk, warm layer and ISAF) obey their 
own equations of continuity and angular momentum con- 
servation, with the appropriate sources terms to account 
for the exchange of matter between the them. The angu- 
lar momentum exchange is given by the mass exchanged 
and the difference in specific angular momentum; viscous 
coupling between the layers is neglected. The continuity 
equation is best expressed as an equation for the local 
mass flux M X (R) (where x stands for either 'ssd', 'warm' 
or 'isaf'). The equation for this is: 

M X (R) = M x {R 0x ) -2k f MR)RdR , (7) 

where Rq x is the inner edge of phase x, which is i?j n for the 
ISAF, i? W arm for the warm layer and i? S sd for the accretion 
disk. The quantity ip x is the mass source term. For the SSD 
one has tp ssd = -0isaf^ ss d - V'ssd^warm, for the warm layer 
one has VWro = Vwi^warm - VWnwisaf and for the ISAF 

One has ?Aisaf = V'warm— >isaf — V'isaf —>ssd • 

The integration constant M x (i?o x ) is zero for the warm 
layer and the SSD, since both have clearly defined inner 
edges due to evaporation. But it is non-zero for the ISAF, 
since this is the layer that in the end carries all the matter 
into the black hole. For the ISAF the value of Misai(Rm) is 
determined as an eigenvalue of the problem: by demanding 
a vanishing mass flux at the outer edge Mi sa f (i? ut) = 0, 
the value of A/j sa f (Ri n ) is indirectly determined. 

In the thin disk approximation the equation for angu- 
lar momentum conservation with a no-friction boundary 
condition at R — Rq x can be cast into the general form: 

" ll M X (R) 



S x ^x 
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In order to solve the S x from this, one must evaluate the 
viscosity coefficient 2/ x , which for a standard a— disk is 
given by 

/,:/; 0) 



a 



i? 3 



GM jim p 



2.4. Numerical approach 



The above equations are solved in a set of nested sub- 
routines and functions. At the deepest level we solve 
the structure of the ISAF for a given evaporation rate 
VWm^isaf (R) and given M sa f {Rin}- This function returns 
the value of the ISAF mass flux at infinity M; sa f(i? — > oo). 
By feeding this function into a root-finding routine (we use 
zbrent from Numerical Recipes) the value of Mi sa f (i?; n ) 
is found for which Mj sa f(i? — > oo) = 0. Note that this 
value is, due to recondensation of the outflowing part of 
the ISAF material, smaller than the total evaporation rate 
2ir J "0warm^isaf (R)RdR. In other words: some material 
goes more than once through the process of evaporation. 
In any case, the above procedure solves the ISAF structure 
for given VWm— isaf (R) 

One step higher in the hierarchy we solve the struc- 
ture of the warm layer, for given values of -R warm and 
i? ss d- First we integrate the equations from i? W arm to i? ss( j 
and determine ^> wa rm— >isaf (R)- We then solve for the struc- 
ture of the ISAF in the way described above. The ISAF 
spreads outward beyond i? ss d and injects its protons into 
the SSD. This proton injection provides the energy needed 
to evaporate the SSD into the warm layer, and therefore 
determines the evaporation rate V'ssd— >warm- This then en- 
ables us to continue the integration of the equations for 
the warm layer from i? ss d out to R — > oo. 

The function that computes the structure of the warm 
layer in the above way returns the value of M warm (R — > 
oo). We feed this function into zbrent so that the value 
of i?warm for which M warm (i? — > oo ) = is automatically 
determined. Note that this is a nested use of zbrent, since 
within the function that determines the structure of the 
warm layer we also call zbrent for determining the struc- 
ture of the ISAF. 
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Fig. 2. The radial three-layered structure of the evaporat- 
ing disk and the hot accretion flows, for the case of low a 
(i.e. a = 0.05) for the warm layer and the ISAF (model 
1). Upper panel shows the surface density of the two com- 
ponents: the cool (SSD) disk, warm surface layer and the 
ISAF. Lower panel shows the local mass flux of the com- 
ponents in units of the Eddington rate (the Eddington 
accretion rate is defined here using an efficiency factor of 
0.1). The horizontal dashed line in the top figure is the 
value of S = 6 g/cm 2 . 



3. Resulting solutions 

Using the above procedure we find stationary solutions for 
most (though not all) parameter combinations. In Figs. [21 
1)1 II we show three examples, all of which have Mbh = 
10 M©, £ = 1 and 77 = 0.9. In model 1 (Fig.EJ) we take a = 
0.05 for all three components, while in model 2 (Fig.|3J the 
ISAF and the warm layer have a ten times higher viscosity: 
a = 0.5. In both these models we take i? ss d = 30 R s . In 
model 3 we take a — 0.5 and i? ss d = 100 R s . As can be 
seen in these figures, the mass flux M for the ISAF is 
positive inward of a certain radius, and negative outward 
of that radius. This negative mass flux is the outflow in 
the ISAF driven by the angular momentum conservation. 
It rapidly goes back to zero toward larger radii due to 
the loss of ISAF ions into the SSD. The mass flux of the 



warm layer is zero at both small and large radii, and has 
a peak value at the inner edge of the SSD (i? ss d)- Beyond 
R ss d the mass of the SSD is evaporated into the warm 
layer, while at R < i? ss d the warm layer evaporates into 
the ISAF. The surface density of the warm layer makes 
a jump at i? ss d, which is due to the temperature jump 
(300 KeV inward of i? ss d and 80 KeV outward of i? ss d)- 
In the case of low viscosity (model 1, Fig. [2J the warm 
layer only narrowly reaches the saturation surface density 
Swarm = S = 6.0, but for the high viscosity case (model 
2, Fig. |2J this saturation value is reached easily. For the 
same high a but with larger i? ss d, the surface density does 
not reach saturation anymore. 

Interestingly, this saturation of S warm causes the mass 
flux in the SSD to acquire an peculiar wiggle in model 2 
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(clearly seen in Fig. [3]): first (at large radii) the matter in 
the SSD is evaporated into the warm layer {dM sst \/dR > 
0), but as one gets closer to the inner radius i? ss d, some 
matter of the warm layer is dumped back onto the SSD 
(dM ss< i/dR < 0). Only just beyond i? ssc j the evaporation 
of the SSD kicks in again and produces the cut-off in M ssc j 
and the maximum in A/ warm at i? ss d- This recondensation 
takes place because if the surface density of the warm 
layer has reached the saturation value So (which is inde- 
pendent of radius), the mass flux in the warm layer should 
go roughly proportional to i? 3 / 2 because both S warm and 
the temperature of the warm layer are then constant with 
R. A decreasing mass flux in the warm layer in the flow 
direction (dM warm /dR < 0) inevitably means reconden- 
sation onto the SSD. This wiggle is not present in the 
case of low viscosity or the case with large i? ss d, since the 
saturation value S is not reached by the warm layer. 

Another interesting difference between the example 
models is that the 300 KeV part of the warm layer (i.e. the 
part inward of i? S sd) is narrower for the case of high vis- 
cosity (model 1) than for the case of low viscosity (model 
2). ' 



3.1. Global trends 

The dependence of some global parameters of the models 
on the input parameter i? ss d are discussed in this subsec- 
tion. In Fig. [5] the width of the 300 KeV part of the warm 
layer (i.e. the part inward of R ss d) is shown as a function 
of i? ss d- It is clear that for low a and for large i? ss d this 
width is the largest. Since this region may be responsible 
for the production of hard X-rays, its surface area is an 
important quantity of the model. Also, if this region is 
wide, then the approximation that the soft photons of the 
SSD do not cool this region (hence allowing it to have 300 
KeV temperature) is more valid than in the very narrow 
case. 

In Fig.|Slthe circulation rate is shown. It tells how often 
the matter is evaporated and recondensed again before 
accreting onto the hole. As one can see, for all values of 
a and i? ss d this rate is always of the order of 3, meaning 
that on average the ISAF matter gets injected into the 
SSD (and re-evaporated from the warm layer) twice before 
finally flowing into the black hole. In other words: 2/3 of 
the hot ISAF matter is used up for the proton injection 
into the SSD. 
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Fig. 5. The width of the 300 KeV region of the warm 
layer as a function of R ss( i- This width is defined as 

-^ssd / ^warm 1 ■ 
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Fig. 6. The circulation rate within the stationary solution. 
3.2. Radiative output 

The three components of the flow each contribute to the 
X-ray spectrum. The cool disk produces a very soft com- 
ponent, the warm layer and the ISAF a hard component 
(peaking above 100 keV). A detailed calculation of the 
overall spectrum produced by the combination of these 
components would be somewhat premature, in view of the 
uncertainties in the model. The relative contributions of 
hard and soft X-rays, however, can be determined sim- 
ply from the energy budgets of these components. Fig. [7| 
shows these in the form of the vertically integrated heating 
rates of the components, per unit surface area. The ISAF 
has the largest amplitude, since the energy release in it 
takes place in the deepest part of the gravitational poten- 
tial. The radiative efficiency of an ISAF is low, however, 
so that the radiative output is actually dominated by the 
other two components of the flow, in particular the warm 
component. The exact value of the radiative efficiency of 
an ISAF depends on details, but typical numbers are of 



Fig. 7. Energetics of the three components of the accre- 
tion flow for model 1 (i.e. a = 0.05 for all components). 
Solid: radiation from the cool disk, dashed: viscous dis- 
sipation in the warm component, dot-dashed: heating of 
the warm layer by ion illumination from the ISAF, dotted: 
viscous dissipation in the ISAF. Since the ISAF is radia- 
tively inefficient, the total radiative output from the flow 
is dominated by the hard radiation from the warm layer. 
The accretion rate for this case is m — 0.06 (in units of 
Eddington) . 



the order e sa ma 2 (e.g. Rees et al. ll982i Esin et al. Il997ft . 
For the case shown in Fig. [7] this yields e « 1.5 x 10~ 4 . 
The contribution of the ISAF is thus small, and most hard 
radiation comes from the warm layer. The total energy 
deposition into the warm layer (viscous heating plus ion 
illumination) is Qtot,warm = 1-2 x 10 35 erg/s for the exam- 
ple of Fig. (model 1). All of this goes into hard radia- 
tion, but a fraction of this radiation is intercepted by the 
cool disk and reprocessed into soft flux. The internal vis- 
cous dissipation in the cool disk is a smaller contribution 
(Qtot,ssd = 4.6 x 10 34 erg/s). In this way, we find, for the 
example of Fig. Fh a rd/F so ft 2.6. The cool disk pro- 
duces a standard multi-color disk spectrum; since its inner 
edge is at some distance from the hole, its contribution to 
the X-ray spectrum is very soft. Most of the energy dissi- 
pated in the warm layer is converted to to hard X-rays by 
Comptonization of soft photons from the cool disk (a frac- 
tion is spent in poducing soft flux by irradiation of the cool 
disk surface). The ISAF also produces a hard spectrum, 
but since it is radiatively inefficient, its contribution to the 
total hard flux is much smaller than its energy budget as 
shown in Fig. 

The conclusion is thus that the model produces a 
predominantly hard spectrum. The soft flux produced is 
mostly due to reprocessing of hard flux by the cool disk 
directly underneath the warm layer, but includes a small 
contribution from internal viscous heating in the disk. The 
hard flux comes mostly from the region where the warm 
layer evaporates into the ISAF, i.e. a ring just inside the 
inner edge of the cool disk. This is consistent with the 
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close connection found observationally between the tim- 
ing properties of the hard flu x and the inner radius of the 
cool disk (e.g. G ilfano v et al. l2OO0t Churazov et al. l200ll 
Revnivtsev et al l200ll) . 

3.3. Stability 

3.3.1. Stability of the vertical structure 

If we envisage the radial structure of the cool disk as fixed, 
one can study whether the structures of the warm and 
hot (ISAF) layers is stable to small perturbations. This 
question was answered in part by Deufel et al. l)2002l) . In 
these calculations, the vertical structure of a cool disk plus 
warm layer under ion illumination was determined in a 
time dependent way, but without the radial structure that 
we have included here. In these calculations, the results 
quickly converge to a stable static equilibrium, as long as 
the cool disk is present. The adjustment takes place on the 
thermal time scale of the warm layer. This stability can be 
understood as a con sequence of the overall energy balance 
(Haardt et al. llflfl^ between comptonization in the warm 
layer and soft photon production in the cool disk. Inside 
the inner edge of the cool disk, however, the warm layer 
becomes unstable due to the lack of soft photon flux, and 
transforms into the evaporating flow studied here. 
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Fig. 8. The inner radius of the cool disk i? ss( j and the in- 
ner radius of the warm phase i? war m as a function of global 
accretion rate M in units of the Eddington accretion rate 
Me (defined using a 10% efficiency factor). The grey area 
shows therefore the extend of the 300 KeV evaporating in- 
ner part of the warm layer inward of the cool disk. Shown 
are the results for models 1 (a = 0.05) and 2 (a = 0.5). 
Note that in the modeling procedure the i? ss d is the model 
parameter and M the result. But since in nature it is the 
opposite, we plot M on the x-axis here. 



3.3.2. Stability of the radial structure 

In the present calculations, additional questions of stabil- 
ity arise related to the radial structure of the accretion 
flows. The first question we posed is if the combination 
of ISAF, warm layer and cool disk is stable if the posi- 
tion of the inner edge of the cool disk is kept fixed. In 
other words, if the evaporation into the ISAF, its spread- 
ing over the cool disk and condensation onto it converges 
to a stable stationary state. Preliminary time-dependent 
simulations of the warm disk, for a fixed accretion disk, 
and for a stationary ISAF computed at every time-step, 
show that the warm layer seems to be stable, and will find 
the stationary solution quickly. But full time dependent 
simulations of both hot phases (the warm disk and the 
ISAF) have not been carried out to confirm this. 

A second question is whether the position of the inner 
edge of the cool disk is stable: on time scales longer than 
the adjustment time of the warm layer, the balance be- 
tween evaporation and condensation may lead to a drift 
of the position of the inner edge of the cool disk. This ques- 
tion can be answered already without a full time depen- 
dent study by considering the relation between the radius 
of the inner edge of the accretion disk and the total mass 
accretion rate. This is shown in Fig. |H1 It shows that for 
increasing accretion rate the inner radius of the cool disk 
(as well as the inner radius of the warm layer) increases as 
well. The reason is that in our model the warm flow near 
i? ss d must carry the entire mass accretion rate M while 
it's surface density is limited to at most So. For higher M 
it therefore requires a larger circumference to transfer all 



the matter with the same surface density, hence the larger 

-Rssd- 

This relation is opposite to what is found from 
the coronal evaporation mechanism of Meyer & Meyer- 
Hofmeister (1994), where the inner radius of the cool disk 
decreases for increasing global accretion rate. This rela- 
tion caused the Meyer & Meyer-Hofmeister evaporation 
mechanism to be self-regulating: if the i? ssc i is too large, 
it will decrease until a balance is reached between the 
total mass accretion rate in the disk and the evapora- 
tion rate. Such a self-regulatory mechanism is not present 
in ion-bombardment evaporation process. A small devia- 
tion from the stationary solution may let the i? S sd move 
e.g. slightly inward, which decreases the allowed global ac- 
cretion rate for a stationary solution. Since the accretion 
rate is fixed from the outside and cannot be changed, it 
will exceed this new value and therefore decrease i? S sd even 
further. We therefore expect that a small deviation of R ss d 
from the equilibrium value would amplify and cause the 
inner edge of the SSD to move inwards or outwards in a 
run-away fashion. The time scale to be expected for this 
instability is the viscous time scale of the SSD. 

The question thus arises whether the evaporating so- 
lutions found here can survive in practice; in any case it is 
clear that some additional piece of physics will be needed. 
This need not be seen as a debilitating disadvantage, how- 
ever, since the variability of the hard state shows that the 
accretion flow is in fact violently unstable. If the evapora- 
tion mechanism had predicted a stable stationary flow, it 
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would have needed additional physics as well, in this case 
something to explain the observed instability of the flow. 

4. Discussion and conclusions 

In this paper we have presented the structure of the 3- 
layere d accretion disk model envisaged by Spruit & Deufel 
( 2002). This models is based on the ion-illumination model 
for X-ray binaries in which a hot ion support ed flo w coex- 
ists with a cool accretion disk (Deufel et al. 120021 : Spruit 
and Deufel l2002h . This coexistence leads to a three- layered 
structure: between the cool disk and the ion supported 
flow a 'warm' layer is formed by ions penetrating the 
disk from the ion supported flow. The ion and electron 
and electron temperatures are in equilibrium in this warm 
layer (in contrast with ion supported flow, where the ion 
temperature is much higher than the lectron temperatue) . 
Its temperature, in the 100 keV range, is governed by the 
equilibrium between ion heating and Compton cooling by 
soft photons from the cool disk. This equilibrium is quite 
stable, in contrast with a single optically thin accretion 
flows in this temperature range, which are strongly ther- 
mally unstable (the well known instability of the 'S LE' 
branch of accretion soluti ons, e .g. Chen et al. [l995). In 
the paper of Deufel et al. (2002) the local vertical struc- 
ture (for a given distance from the hole) of the accretion 
flow was calculated. The calculations presented here are 
complementary to these. They focus instead on the radial 
structure of the accretion flow, while the vertical structure 
is approximated by a 3-layer approximation. The exchange 
of energy and angular momentum between these layers is 
taken into account, as well as the processes of evaporation 
and condensation between the ion supported flow and the 
warm layer, and between the warm layer and the cool disk. 

Wc find that stationary solutions exist for a wide range 
of parameters, and that these flow geometries exhibit cir- 
culation patterns in which matter is evaporated and re- 
condensed several times before vanishing into the black 
hole. This finding is of special significance for the Li- 
enhancements observe d in the secondaries of black hole 
binaries (Martin et al. ll994h . which can be explained as a 
result of Li-production in the accretion disk by ion illumi- 
nation (spallation of CNO nuclei and He+He fusion), pro- 
vided that in addition to the accretion flow there is an out- 
flow from the inner regions of the accretion disk that can 
carry the Li to the secondary. This is an attractive possi- 
bility since the hard spectral states that we associate here 
with the ion-illumination process is known to be closely 
related to the presence of jets in in black hole binaries. 
The 'recirculation' of matter through the ion-bombarded 
warm layer increases the Li-production significantly. 

We find that the energy radiated by the ion- 
illuminated, evaporating flow is likely to be dominated 
by a hard X-ray component produced in the evaporat- 
ing region itself. A soft component is also produced by 
reprocessing in the cool disk, but it is energetically less 
important than the hard component. This provides a sat- 
isfactory solution to a long-standing problem in models 



where the hard and soft fluxes are coupled through repro- 
cessing. The generic model for this coupling (Haardt and 
Maraschi) has the property of predicting similar energy 
fluxes in the hard and soft components. In order to make 
the hard component dominate as observed, modifications 
of the geometry are needed to reduce the interaction be- 
tween the soft photons and the Comptonizing region. In 
the evaporating flows presented here, this occurs naturally 
because most of the hard flux is produced in a 'soft photon 
staved' region just inside the inner edge of the cool disk. 

We also show that the globally stationary solutions are 
not stable on long time scales. Eventually the inner rim of 
the evaporating SSD will move either outward or inward 
in a run-away fashion. However, if the inner rim of the 
SSD is kept fixed, the structures of the warm intermediate 
layer and the hot ISAF are presumably stable. Hence, this 
evaporation mechanism is feasible, but it will result in 
time-dependent phenomena instead of a stationary state. 
How these time-dependent solutions behave remains the 
topic of a future detailed time-dependent study of ion- 
illuminated evaporating disks. 
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